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Program Summary: 

• Program name: KROWIG. 

• Version: 1.0 (June 1991). 

• Author: Thorsten Ohl, <ohl(9crunch. ikp .physik. th-darmstadt . de> 

• Long write-up: This document. 

• Programming language: FORTRAN-77. 

• Computer/operating system: Any with a FORTRAN-77 environment. 

• Number of program lines: ^ 3000 

• Other programs used: KRONOS0 (Version 1.2 or later), HERWIGpOj] (Ver- 



sion 5.4), PAKPDFp 



• Input files needed: None. 

• Initial parton shower: HERWIG. 

• Hard subprocesses generated: Neutral current deep inelastic scattering. 

• Final parton shower: HERWIG. 

• Fragmentation model: Cluster (HERWIG). 

• Initial QED radiation: Leptonic: leading logarithmic approximation, 
summed to all orders (KRONOS); quarkonic: N/A. 

• Final QED radiation: Leptonic: leading logarithmic approximation, 
summed to all orders (KRONOS); quarkonic: included in the HERWIG par- 
ton shower (optional). 
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1 Introduction 



The new ep collider HERA will explore a new kinematical range for deep inelastic 
lepton-hadron scattering (DIS) Due to the smaller mass of the electron and 
the higher available energy, radiative corrections will be much larger than at the 
previous /ip experiments. Inclusively, the leading contribution in each order of 
perturbation theory can be estimated as 

(1) 

and the exclusive corrections can surpass 100% in some regions of phase space 0. 
A quantitative understanding of the radiative corrections is therefore mandatory 
for the physical interpretation of HERA data. The electroweak corrections to 
deep inelastic ep scattering have been calculated completely to one loop order 

^] and the leading logarithms of the next order have been given 0. 

However, for detailed experimental studies, the implementation of these cal- 
culations in a Monte Carlo event generator is indispensable. The order uqed 
corrections have been implemented in the HERACLES event generator 0, which 
has been interfaced in DJANGO0 to the LEPTO QCD Monte Carlo [g. At the 
parton level, the leading higher order QED corrections have been implemented 
in the KRONOS Monte Carlo [|]. 

KROWIG extends the hadronically inclusive predictions of KRONOS to exclusive 
hadronic final states. In addition, KROWIG allows to incorporate radiative cor- 
rections into the predictions of the HERWIG model for QCD parton showers and 
cluster fragmentation [1^. Thus KROWIG is unique in two different aspects of 



HERA physics: it is the first Monte Carlo Event generator implementing QED 
radiative corrections to DIS in the HERWIG framework and it is also the first Monte 
Carlo event generator implementing higher order QED radiative corrections for 
hadronically exclusive final states. Although multiphoton events will not play 



an important role during the first years of experimentation at HERA |Tl|, the 
included soft photon exponentiation can have numerically appreciable effects. 

After elaborating on the physics of radiative corrections to exclusive hadronic 
final states in section we shall discuss their implementation in KROWIG in 
section p.2| . The various parameters controlling KROWIG are discussed in section |[ 



Section ^ is devoted to the FORTRAN-77 interface for application programs. Before 
concluding, we list in section ^ the limitations of KROWIG version 1.0. Technical 
details of the installation of KROWIG and a sample application are relegated in the 
appendix. 

In order to avoid unnecessary duplication, this manual frequently refers to the 
manuals for KRONOS |12| and HERWIG [jl3[ and assumes familiarity with the two 



programs. 
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2 Implementation of KROWIG 1 . 



2.1 Physics issues 

Conceptionally, the interfacing of leptonic QED radiative corrections and QCD 
Monte Carlos is simple, provided one stays in the framework of the leading log- 
arithmic approximation. In this case the radiative corrections factorize and can 
be presented to the QCD Monte Carlo as just another hard subprocess. 

In a nutshell, the reasoning is that the hadronically inclusive cross section 
including radiative corrections is given by the parton model, while QCD parton 
showers and the subsequent fragmentation must happen with unit probability. 

Factorization is in our approach equivalent to staying within the usual one- 
photon-approximation, on which the QCD description of DIS in terms of the 
improved parton model or operator product expansion is built. In particular 
this means that we will neglect the non-leading box diagrams in which a second 
photon (or Z°) is exchanged between the leptonic and hadronic subsystems. 

Using this approximation, the hadronically inclusive cross section can be writ- 
ten as 

"^'"^ i^,y,z) = ^^E^7.(x',g'^)(A. + (i-y05.) (2) 



dxdydz x'^y''^s''^ 



i=q,q 



where 



X = — ■ -, y = , s = zs, Q = zQ (3) 

y + z — 1 z 

In these formulae x and y are the (leptonic) Bjorken variables and z is the energy 
fraction of the incoming electron after initial state radiation. Final state radiation 
has been ignored here for simplicity, but is included in KROWIG. Ai and Bi denote 
the appropriate combination of electroweak couplings and propagator factors for 
the quark i. D and the fi are solutions of the familiar evolution equations 



dlnQ"^ J rj 3 271 \1 — rj 



dD{^,Q^) _ rdriaQED l + ri 



2' 



m/v,Q') (4) 



DU/V,Q') (5) 



(6) 

which are summing the leading logarithms (P and their QCD counterparts, re- 
spectively. 

The Monte Carlo implementation of (^) can now be decomposed into five 
parts: 
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• Generation of {x,y,z) triples according to the distribution (g). 

• Generation of photons corresponding to the chosen {zjQ"^). 

• Generation of final state photons. 

• Generation of QCD parton showers corresponding to the chosen [x'^Q'"^). 

• Hadronization of the generated partons. 

It is now obvious that these problems can be solved by almost independent pro- 
grams. Because the hadronically inclusive cross section is known, the first three 
points can be handled by KRONOS without any intervention from HERWIGQ Sim- 
ilarly, HERWIG needs no information from KRONOS besides the (x', Q'"^) pair to 
reconstruct an appropriate parton shower. Furthermore, as already stated above, 
in the parton model hadronization happens with unit probability, once the hard 
subprocess has been generated. 

From this discussion a simple physical picture emerges: the Monte Carlo event 
generator for radiative corrections at the parton level can be used unchanged, be- 
cause it only depends on the hadronically inclusive cross sections. The algorithms 
implementing this step in KRONOS have been described in . On the other hand, 
the QCD Monte Carlo can completely ignore the leptonic radiative corrections, 
provided it uses the electron after initial state radiation as input. Phrased differ- 
ently, the leptonic radiative corrections affect the QCD Monte Carlo only through 
broadening the sharp electron energy spectrum into an energy distribution. 

2.2 Programming issues 

In order to exploit the simple physical picture of initial state radiation as an 
effective beam energy spread, several requirements must be met. 

First of all, the called QCD Monte Carlo event generator must be very dis- 
ciplined in its use of global variables. Any assumption based on a fixed beam 
energy will defeat a simple interface. Typically such assumptions would be used 
at initialization time to set up internal tables, whose contents will be invalid for 
radiative events. 

Fortunately, HERWIG is a well-behaved program in these respects. Since ini- 
tial state radiation can only reduce the effective beam energy, HERWIG's table of 
Sudakov formfactors is not invalidated by radiative events. We can therefore use 
the following simple multiply layered approach, which is shown graphically in 
figure |1]. 

-'^ There is one subtlety, however: when HERWIG reconstructs the initial state parton shower 
from the chosen (x', Q'^) pair, it might reject it due to mismatches between the ordering in 
the evolution equations of the employed parton distribution functions and the angular ordering 
in HERWIG's evolution. This is taken care of in KROWIG by rejecting the event. 
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KRONOS driver program and parameter management 



krowigO HWIGPR) 



HWEPROO 



HWHDISO 



kronos 



Figure 1: Top level structure of KROWIG. 



The main routine is derived from the command interpreter and parameter 
management routines of KRONOS, which have been described in in detail. The 
only difference is the larger number of parameters which are needed to con- 
trol HERWIG| (cf. table Instead of calling the parton level kronos () event 
generation routine, the main routine now calls the krowigO routine which has 
been adapted from HERWIG's main event generation loop. In the innermost layer, 
HERWIG's parton level generator is emulated by a wrapper around the original 
kronos parton level generator. 

This simple structure is straightforwardly implemented because both pro- 
grams use the standard /hepevt/ event record [jl4|, which allows a transparent 
flow of event information between Monte Carlo components. Otherwise more 
complicated translation routines would have been necessary. 



3 Parameters 



3.1 KRONOS Parameters 



A detailed discussion of the KRONOS specific parameters can be found in & 12 



and will not be repeated here. The parton distribution functions are taken from 



^It is important to note that this implies that the parameter common blocks /krpcom/ of 
KRONOS and KROWIG are not compatible. 
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Variable name 


semantics 


Default value 


ahpla 


'^/otQED 


137.0359895 


lambd 


^QCD 


0.18 GeV 


massle 


me± 


0.51099906 • 10-3 GeV 


masslp 


rup 


0.93827231 GeV 


masslz 


Mzo 


91.17 GeV 


masslw 


Mw± 


80.0 GeV 


sin2w 


sin^ 9w 


0.23 GeV 


elecen 


e~ energy 


30 GeV 


proten 


p energy 


820 GeV 


nevent 


Number of events 


1000 


bstyle 


"branching style"0 


4 


epsiln 


Internal infrared cutoff 


10-^ 


cutxmi 


Lower cut in x 


0.001 


cutxma 


Upper cut in x 


0.5 


cutymi 


Lower cut in y 


0.1 


cutyma 


Upper cut in y 


0.9 


cutqmi 


Lower cut on hadronic 


100 GeV2 


cutqma 


Upper cut on hadronic 


s 


rseed 


Random number seed 


54217137 


errmax 


maximum error count 


100 


verbos 


verbosity 





run id 


run identification 




stdin 


standard input 


5 


stdout 


standard output 


6 


stderr 


standard error 


6 



Tor reasons of backward compatibility, this parameter has slightly confusing semantics 
associated with it. It should be noted however, that the value of bstyle is not at the discretion 
of the user. The only supported values are 

4 Full radiative corrections switched on (default). 

No radiative corrections. Born cross section. 

All other values will be reset to the default value. 

Table 1: KRONOS parameters 
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the PAKPDF library |]T5[ and can be set with the usual KRONOS commands, which 
are also explained in detail in P, |12|. Note however, that the optional use of 
phenomenological parametrizations of electroproduction structure functions Fi 
is not supported, because KRONOS needs parton distributions for generating the 
struck quark. 



3.2 HERWIG Parameters 

The variables in table ^ fall into three categories 

• herwig and soe have no analogue in HERWIG, they control whether and how 
HERWIG is called to generate hadronic final states from KRONOS' partons. 

• some HERWIG variables have denote the same physical parameter as a KRONOS 
variable; here the former is set from the latter. 

• finally there are plenty of variables in HERWIG that have no analogue in 
KRONOS: the name of these variables is derived from their name in HERWIG 
by prepending 'hw' and stripping the fifth and sixth character (this does 
not lead to ambiguities for the variables considered). 

HERWIG variables which do not affect DIS have not been made available in the 
krdcmdO interface. 

Note that the variable VPCUT relates to photons radiated from the hadronic 
subsystem and does not affect the leptonic side, which is handled by KRONOS. 



3.3 Tuning of parameters 

In the Monte Carlo generator working group of the 1991 HERA workshop 
HERWIG version 5.3 has been tuned to EMC data, using the KMRS BO parton 
distributions. The following two sets of changes with respect to the HERWIG 
version 5.3 default values have been proposed: 

1. No soft underlying event: 



# herawglO-setl .krowig 

set sue false 

set hwclma 2.00 

set hwqspa 2.00 

set hwpspl . 50 

set hwptrm . 70 
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variable 


HERWIG 


semantics 


Default value 


herwig 


N/A 


call KROWIGQ 


true 


sue 


N/AQ 


generate S.U.E. 


false 


elecen 


PBEAMl 


electron beam energy 


30 GeV 


proten 


PBEAM2 


proton beam energy 


820 GeV 


nevent 


MAXEV 


number of events 


1000 


lambd 


QCDLAM 




0.180 GeV 


hwvqcu 


VQCUT 


quark virtuality cutoff 


0.48 GeV 


hwvgcu 


VGCUT 


gluon virtuality cutoff 


0.10 GeV 


hwvpcu 


VPCUT 


photon virtuality cutoff 


10^0 GeVg 


hwclma 


CLMAX 


maximum cluster mass 


3.35 GeV 


hwpspl 


PSPLT 


cluster mass distribution^ 


1.00 


hwqdiq 


QDIQK 


max. scale for glue —>■ diquark 


0.00 GeVg 


hwpdiq 


PDIQK 


glue diquark rate 


5.00 


hwqspa 


QSPAC 


spacelike evolution cutoff 


2.50 GeV 


hwptrm 


PTRMS 


intrinsic transverse momentum 


0.00 GeV 


hwenso 


ENSOF 


multiplicity of S.U.E.Q 


1.00 


hwipri 


IPRINT 


HERWIG's verbosity^ 





hwmaxp 


MAXPR 


# of events printed 





hwmaxe 


MAXER 


maximum # of errors 


100 


hwlwev 


LWEVT 


unit for writing output events^ 


oE 


hwlrsu 


LRSUD 


unit for restoring formfactors 





hwlwsu 


LWSUD 


unit for saving formfactors 


77 


hwazso 


AZSOFT 


soft gluon correlations 


true 


hwazsp 


AZSPIN 


gluon spin correlations 


true 


hwncol 


NCOLO 


# of colors 


30 


hwnf la 


NFLAV 


^ of active flavours 
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"Iff this variable is false, then KROWIG behaves identical to KRONOS. 

''This is implemented by adding 10000 to the HERWIG variable IPROC iff sue is false. 

'^i.e. photons in hadronic final state parton showers are suppressed by default. 

''See for a detailed description of this parameter. 

'^i.e. no diquarks by default. 

■''See [13 1 for a detailed description of this parameter 

^See |1S| for a description of this parameter. The default value shuts up HERWIG almost 
completely. 

''This option should not be used because radiated photons might be missing and HERWIG's 
internal coordinate system is different from KROWIG's (KROWIG has the incoming proton in the 
+z direction, whereas HERWIG assumes that the incoming proton goes in the —z direction). 

*i.e. don't write events to disk. 

^ You probably don't want to change this! 



Table 2: HERWIG parameters accessible through krdcmd ( ' set 



') 
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* krowigappl . f 

call krdcmd ( ' set maxpr 1 ' ) instruct HERWIG to dump the first event 

call krdcmd ( ' init ' ) initialize the generator 

call krdcmd ('generate 10000') generate 10000 events 

call krdcmd ( ' close ' ) cleanup 



Figure 2: FORTRAN-77 interface 

In addition, it has been proposed to set the parameter BTCLM in the subrou- 
tine HWCCUTO to 3.00, while increasing CLMAX to 3.00. However KROWIG does 
not support the modification of BTCLM because it is not a tunable HERWIG 
parameter and changing it requires modification of the HERWIG sources. 
This should be done by the user at his own risk. If 

2. With soft underlying event: 

# herawgl0-set2.krowig 
set sue true 
set hwclma 2.00 
set hwptrm 0.40 
set hwenso . 60 



This set gives a slightly worse agreement with EMC data. 

It should be noted however, that the default parameter values have been 
changed from HERWIG 5.3 to HERWIG 5.4. Therefore a retuning might be necessary. 

4 FORTRAN-77 Interface 

The application program interface of KROWIG is almost identical to KRONOS' Q. On 
the top level, where KROWIG is controlled by character strings passed to krdcmd () 
(see figure II), the interface is identical, except for the possibility to change HERWIG 
parameters. These parameters are collected in table 0. The generate command 
will execute a loop similar to the one shown in figure |^. The generated event will 
be passed to the application program in the standard /hepevt/ common block 
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* krdcmd.f 

subroutine krdcmd (cmdlin) 
character* (*) cmdlin 



else if (cmdlin. eq. 'gen 
do 10 n = 1 , nevent 

call krowig (1) generate an event 

call hepawk ('scan') analyze the event 

10 continue 
else 

end 



Figure 3: Event generation loop 



At the lower level function call interface there are also almost no changes 
from KRONOS. The call to the kronosO subroutine is to be replaced by a call to 
the wrapper routine krowig () which calls kronosO and the appropriate HERWIG 
routines. The single integer parameter is interpreted as follows: 

0: initialize the generator and write an initialization record to /hepevt/. 

1: generate an event and store it in /hepevt/. 

2: perform final calculations and write the results to /hepevt/. 

As in any interface of independly developed programs, there remain some 
rough edges. While these rough edges could be avoided in areas where accepted 
standards exist (e.g. the passing of event information via /hepevt/ fl^ is by 



now well established), the control of input parameters still leaves something to 
be desired: the fine grained control over reinitializations in KRONOS (see for 
a description) can not be achieved with HERWIG, unless we allow to make some 
(simple) modifications to the HERWIG sources. Because it is desirable to keep these 
untouched^, we let modifications of all HERWIG parameter cause a reinitialization, 
even if not strictly necessary. 



5 Limitations of KROWIG 1.0 



•^This should allow KROWIG to work with future HERWIG versions unmodified. 
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5.1 Compton channel 

KROWIG is not applicable in the kinematical region corresponding to the so-called 
Compton events: 

(Pe+P7)'»QL*.^0 (7) 

In this region the quark parton model is not apphcable and therefore HERWIG does 
not provide a realistic model for the generation of hadronic final states. For stud- 
ies of the leptonic and photonic final states of these events it is recommended to 
use the stand-alone KRONOS with the structure function {not parton distributions) 
option. 

Future revisions of KROWIG might overcome this restriction by exploiting a 
simple model (maybe along the lines of the HERWIG remnant fragmentation - 
a.k.a. "soft underlying event"). 

5.2 0{as) matrix elements 

KROWIG can not yet be used with the 0{as) matrix element option of recent 
HERWIG versions. This is not a fundamental limitation and could be overcome 
in future versions of KROWIG, after the 0{as) QCD inclusive cross sections have 
been implemented in KRONOS. 

5.3 Charged currents 

No charged current events will be generated. This restriction will be lifted (with- 
out any changes to KROWIG), once the charged current subsystem of KRONOS has 
been released. 

6 Conclusions 

We have presented version 1.0 of the Monte Carlo event generator KROWIG for 
deep inelastic scattering at HERA energies. KROWIG acts as an interface for the 
QED generator KRONOS and the QCD generator HERWIG. It is suitable for the 
study of hadronic final states, taking into account the bulk of the electro-weak 
radiative corrections. 

KROWIG is complementary to the DJANGO Monte Carlo, because it includes 
higher order QED corrections and uses a different QCD Monte Carlo. 

Availability 

The latest release of KROWIG is available from the author upon request. In addi- 
tion, KROWIG is available by anonymous ftp from freehep.scri.fsu.edu in the 
directory f reehep/event_generators/krowig. It is nevertheless recommended 
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to notify the author, in order to be informed of future bug fixes and enhance- 
ments. 
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A Revision History 

Version 1.0, June 1992 

First official release. 

B Installation 



KROWIG is distributed in PATCHY ^ CARDS format. The installion scripts (a 
Makefile for UNIX, a BUILD.COM DCL ffie for VMS, and a JCL/370 deck for 
MVS) assume the availability of the KRONOS CARDS file (version 1.2 or later) and 
of the HERWIG0 (version 5.4) sources in CARDS format. Furthermore PAKPDF Jl5| 
and hepawk [[1^, |19| are expected as linkable libraries. The latter can be replaced 
by any other event analysis subroutine. 

If PATCHY is not available, a FORTRAN-77 version of the source can be pro- 
vided by the author. 

B.l Building HERWIG 

In a ffist step, HERWIG should be compiled without the functions that will be 
overwritten by KROWIG: hwhdis, pdfset, and structf. This can be done by 
removing them from the FORTRAN-77 distribution or by using the CARDS file 
with the CRADLE 

+EXE. 

+USE,IBM,IBMMVS. 
+USE,*HERWIG. 

+USE , P=HERWIG , D=HWHDIS , T=INHIBIT . 
+USE , P=HERWIG , D=STRUCTF , T=INHIBIT . 
+USE , P=HERWIG , D=PDFSET , T=INHIBIT . 
+PAM,11,T=CARDS. 



"^Usually available from CBHEP: :DISK$USER1: [THEORY. HERWIG. PATCHY] HERWIG54. CAR via 
DECNET. 
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+QUIT. 



If the standard HERWIG distribution runs on the target system, no changes should 
be necessary for using it in KROWIG. 

B.2 Building KROWIG 

After adjusting filenames in the installation scripts, these scripts should build 
KROWIG without further user intervention. If HERWIG is not available in CARDS 
format, the HERWIG common blocks have be imported manually into the routines 
requesting them with +CDE , HERCDES. Note that it is not possible to link KROWIG 
with the standard KRONOS object module because the parameter common block 
/krpcom/ has to be enlarged and the call to the kronosO subroutine has to be 
replaced by the krowigO wrapper. 

C Common Blocks and Subroutines 

To avoid possible name clashes with other packages, all external symbols exported 
by KRONOS begin with the two letters KR. This convention is also followed by 
KROWIG, except where routines from other programs are overwritten. 

• Common Blocks: 

The following common block is used by KROWIG in addition to KRONOS' and 
HERWIG's common blocks. 

— /krcevt/: small event record in the style of /hepevt/ for hiding the 
photons from HERWIG. 

— /krcwig/: a collection of variables used for non-local communication 
between KRONOS, HERWIG, and KROWIG. 

• Primary entry point: 

— krowig: main entry point, replacing kronos (which it calls itself). 

• Utility routines: 

— kr2hwl: glue routine, editing KRONOS' output in /hepevt/ to match 
HERWIG's expectations. 

— kr2hw2: glue routine, fixing /hepevt/ after HERWIG has processed it. 

— kr2ent: utility routine, used in kr2hwl for filling /hepevt/. 

• Overwritten entry point: 
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— hwhdis: overwritten standard HERWIG deep inelastic subprocess. 

— hwaend: overwritten HERWIG user exit. 

— pdfset, structf: emulation of PDFLIB 
by PAKPDF [|l|]. 

D Example 

As an example application of KROWIG, we demostrate how the generate a rapidity 
distribution and the so-called "seagull-plot" . 

D.l sample .krowig 

Here is a simple KROWIG command file, setting up parameters and generating 1000 
events. 

# sample . krowig — Example KROWIG 1.0 command file 
# 

# nominal HERA energies 
set elecen 30.0 

set proten 820.0 
# 

# no. of events 
set nevent 1000 
# 

# of course we want QCD included 
set herwig true 

set hwmaxe 100 
set hwmaxp 
set hwlrsu 
set hwlwsu 
# 

# Morf in/Tung 'B2' 
struct mt 4 

# 

# kinematical cuts 
set cutqmi 20 . 
set cutymi 0.015 
set cutyma . 99 
set cutxmi 0.7e-4 
set cutxma . 99 

# 

# go! 
init 

generate 



20 routines (used in HERWIG) 
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# 

# we ' re done . 
exit 



D.2 sample .hepawk 

This is a simple hepawk script, selecting events that satisfy semi-realistic ac- 
ceptance cuts and filling histograms for charged rapidity distributions and the 
famous "seagull plot"0. 

# sample .hepawk — sample HEPAWK analyzer for KROWIG vl.O 

BEGIN 
{ 

printf ("\nWelcome to the KRDWIG test:\n"); 

printf ( " ****************************\n\n" ) ; 

printf ("Monte Carlo Version: °/,s\n" , REV); 
printf (" Run: 7,d\n" , RUN); 

printf (" Date: y.s\n\n" , DATE); 

x_min = 0.0001; 
y_min = 0.02; 
Q2_min = 25.0; 
W2_min = 0.0; 

theta_min_had = theta_min_em = 6.0 / DEG; 
theta_max_had = theta_max_em = 174.0 / DEG; 

h_rap_cms = bookl (0, "rapidity (cms)", 50, -10, 10); 
h_sg = bookl (0, "Seagull", 50, -1, 1); 
h_xf = bookl (0, "xf", 50, -1, 1); 

incut =0; # initialize counter 

> 



S = (@Bl:p + @B2:p)"2; 
f ound_electron = 0; 

for (@e in ELECTRONS) # Collect the outgoing electron. 
{ 

$q = @Bl:p - @e:p; 
Q2 = - $q'-2; 

^Note that this hepawk script requires version 1.2 or later of hepawk ||l^, because the Lorentz 
boost operator / 1 is not available in version 1.0, which has been described in p8[. 
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X = Q2/(2*$q*@B2:p) ; 

y = Q2/(x*S) ; 

if (q2_min <= Q2 && x_iiiin <= x && y_min <= y 

&& theta_iiiin_em <= angle (@e:p, @B2:p) <= theta_iiiax_eiii) 

{ 

f ound_electron = 1; 
break; 

} 

} 

if ( ! f omid_electron) 
next ; 

} 



# Step 1: find hadronic center of mass system: 

$p_hadr = $NULL; # Collect (charged) hadronic momentum 

for (@p in HADRONS) 
if (charge (Op: id)) 

if (theta_min_had <= angle (@p:p, ©B2:p) <= theta_max_had) 
$p_hadr += @p:p; 



W2 = $p_hadr"2; 
if (W2 < W2_min) 
next ; 

incut ++; # We've passed all inclusice cuts, count this event! 

if ($p_hadr:0 ==0) # no hadrons in the cuts, 

next ; # skip this event 

W = sqrt (W2); 



# Step 2: boost to hadronic CMS 
$beta = $p_hadr / $p_hadr:0; 

$q_phot = ($p_hadr - @B2:p) /| $beta; # virt. photon in hadr. CMS 

$q_phot -= $E0 * $q_phot:0; # - make it a 3-vector 

$q_phot /= sqrt (- $q_phot"2) ; # - and normalize 

for (@p in HADRONS) 
if (charge (@p:id)) 

if (theta_min_had <= angle (@p:p, ©B2:p) <= theta_max_had) 
{ 

$p = (§p:p /I $beta; # momentum in hadr. CMS 
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fill (h_rap_cms, rap ($p)); 



xf = - 2 * $p * $q_phot / W; # project on virt . photon 
# extract transversal component and square it 
$p:0 = 0; 

pt2 = - ($p + $q_phot * ($q_phot * $p))"2; 
fill (h_sg, xf, 0, pt2) ; 
fill (h_xf, xf); 

> 

} 

END 
i 

# Dump some numbers 
printf ("\nRESULTS:\n") ; 
printf ("********\n\n") ; 

printf (" raw events: °/od, raw cross section: °/„g mb\n" , 

NEVENT , XSECT) ; 

printf (" error: 7,gn\n" , (ERROR/XSECT) * 100.0); 

printf ("events after cuts: °/od, cross section after cuts: y„g mb\n" , 

incut, (incut/NEVENT) * XSECT); 
printf ( " \nHISTOGRAMS : \n" ) ; 
printf ("***********\n\n") ; 

scale (1/NEVENT) ; 

# finish the seagull plot: 
arith (h_sg, h_sg, "/", h_xf ) ; 
delete (h_xf ) ; 

plot ; # plot the histograms 
printf ( " \ndone . \n" ) ; 

> 

D.3 sample . output 

The following output should result from the input files above, modulo small 
roundoff errors. Note that the warning messages from HERWIG correspond to the 
mismatches in the QCD evolution mentioned in the footnote on page ^ and are 
harmless. 

krdcmd: message: starting KRONOS, Version 1.02/00, (build 920626/0055) 
hepawk: message: starting HEPAWK, Version 1.02/00, (build 920623/0117) 
krowig: message: starting KROWIG, Version 0.99/00, (build 920626/0055) 
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Welcome to the KRDWIG test: 



Monte Carlo Version: vOO.99 (Jun 26 00:00:00 1992) 
Run: 1035996352 
Date: Jun 26 00:56:00 1992 



HERWIG 5.4 



JANUARY 1992 



PARTICLE TYPE 21=PI0 SET STABLE 



INPUT EVT WEIGHT 
INPUT MAX WEIGHT 



, lOOOE+01 
, lOOOE+01 



HWWARN CALLED FROM SUBPROGRAM HWSBRN: CODE = 105 
EVENT 95: SEEDS = 382793286 & 1741717887 

EVENT KILLED. EXECUTION CONTINUES 



WEIGHT = 



lOOOE+01 



HWWARN CALLED FROM SUBPROGRAM HWSBRN: CODE =101 
EVENT 228: SEEDS = 1689384423 & 1215685145 
EVENT KILLED. EXECUTION CONTINUES 



WEIGHT = 



lOOOE+01 



HWWARN CALLED FROM SUBPROGRAM HWBGEN: CODE = 100 
EVENT 310: SEEDS = 929482140 & 326931761 
EVENT KILLED. EXECUTION CONTINUES 



WEIGHT = .lOOOE+01 



HWWARN CALLED FROM SUBPROGRAM HWSBRN: CODE = 105 
EVENT 360: SEEDS = 255018931 & 1779155573 
EVENT KILLED. EXECUTION CONTINUES 



WEIGHT = 



lOOOE+01 



HWWARN CALLED FROM SUBPROGRAM HWSBRN: CODE = 105 
EVENT 381: SEEDS = 915725380 & 739707742 
EVENT KILLED. EXECUTION CONTINUES 



WEIGHT = 



lOOOE+01 



HWWARN CALLED FROM SUBPROGRAM HWSBRN: CODE = 105 
EVENT 519: SEEDS = 461210881 & 1073328688 

EVENT KILLED. EXECUTION CONTINUES 



WEIGHT = .lOOOE+01 



HWWARN CALLED FROM SUBPROGRAM HWBGEN: CODE = 100 
EVENT 690: SEEDS = 281115117 & 1525150099 
EVENT KILLED. EXECUTION CONTINUES 



WEIGHT = .lOOOE+01 
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HWWARN CALLED FROM SUBPROGRAM HWSBRN: CODE = 105 

EVENT 805: SEEDS = 485884270 & 1401730461 WEIGHT = .lOOOE+01 
EVENT KILLED. EXECUTION CONTINUES 
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RESULTS : 
******** 



raw events: 
error: 

events after cuts: 



1000, raw cross section: .4587E-04 mb 

.4994 ■/. 

720, cross section after cuts: .3303E-04 mb 



HISTOGRAMS : 



rapidity (cms) 
HBOOK ID = 



DATE 



8 


4 


9 


8 


2 


X6 


g 




8XX 


7 


Q 


XXX 


J 


6 


XXX 


7^ 


_4 


XXX 


7^ 


2 


9XXX 


7 




YYYYC: 


6. 


.8 


YYYYY 


6 


6 


YYYYY 

AAAAA 


6. 


_4 


YYYYY 
AAAAA 


6. 


_2 


YYYYY 

AAAAA 


Q 




YYYYY 

AAAAA 


5. 


.8 


1 AAAAA Vif 


5. 


.6 


AAAAAAA 


5. 


.4 


xxxxxxx 


5. 


.2 


xxxxxxx 


5 




xxxxxxx 


4. 


.8 


xxxxxxx 


4. 


.6 


9XXXXXXX 


4. 


.4 


xxxxxxxx 


4. 


.2 


xxxxxxxxo 


4 




xxxxxxxxx 


3. 


.8 


xxxxxxxxx 


3. 


.6 


xxxxxxxxx 


3. 


.4 


9XXXXXXXXX 


3. 


.2 


XXXXXXXXXX 


3 




xxxxxxxxxx 


2. 


.8 


XXXXXXXXXX 


2. 


.6 


xxxxxxxxxx 


2. 


.4 


xxxxxxxxxx 


2. 


.2 


xxxxxxxxxx 


2 




xxxxxxxxxx 


1. 


.8 


0XXXXXXXXXX4 


1. 


.6 


xxxxxxxxxxxx 


1. 


.4 


xxxxxxxxxxxx 


1. 


.2 


xxxxxxxxxxxx 


1 




xxxxxxxxxxxx 




.8 


6XXXXXXXXXXXX3 




.6 


XXXXXXXXXXXXXX 




.4 


6XXXXXXXXXXXXXX 




.2 


15XXXXXXXXXXXXXXX920 



CHANNELS 10 1 2 3 4 5 

1 12345678901234567890123456789012345678901234567890 

CONTENTS 1 . 134577886541 

*10** 1 00000000000000001376357193197066100000000000000000 

00000000000000031331995979302197941000000000000000 



LOW-EDGE 



10 
1. 




09988877666554443322211 1122233444556667788899 
06284062840628406284062840482604826048260482604826 



ENTRIES 
BIN WID 



13490 
.4000E+00 



ALL CHANNELS = .6745E+01 
MEAN VALUE = .5607E-01 



* UNDERFL 

* R . M . 
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Seagull 



HBOOK 



ID = 



DATE 



4^ 


.6 


4^ 


.5 


4^ 


_4 


4 


3 


4 


2 


4 




4 




3. 


9 


3. 


8 


■a 

o . 


7 


3 


6 


3. 


.5 


3. 


_4 


3 


3 


3. 


2 


3 




3 




2. 


.9 


2 


8 


2 


7 


2. 


.6 


2. 


.5 


2. 


.4 


2. 


.3 


2. 


.2 


2. 


.1 


2 




1. 


.9 


1. 


.8 


1. 


.7 


1. 


.6 


1. 


.5 


1. 


.4 


1. 


.3 


1. 


.2 


1. 


.1 


1 





.8 
.7 
.6 
.5 
.4 
.3 
.2 
.1 




X 
X 
X 
9X 
6XX 
XXX 



7 2 
X X 
X X 
4X X 
XX X 



3 

X8 
XX 

XX 5 



XX X55XX X 
XX XXXXX X71 
XX3XXXXX XXX 
8XXX X5XXXXXXXX9XXX6 
XXXX0XXXXXXXXXXXXXXX7 
2XXXXXXXXXXXXXXXXXXXXX7 3XXXXXXXXXXXXXXXXXXX X 
XXXXXXXXXXXXXXXXXXXXXXX8329XXXXXXXXXXXXXXXXXXXX1X 
XXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 



6 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 
X 

X 
X X 
X X 
X X 
X X 
X X 
X X 
X X 
X X 
X X 
X X 
X X 
X X 
IX X 
XX X 
XX X 
XX X 
29XX X 
X XXXX X 
X IXXXX X 
X8XXXXX X 
9XXXXXXX X 
7 XXXXXXXX X 
2 X8XXXXXXXX X 
X XXXXXXXXXXIX 
X9XXXXXXXXXXXX 
049XXXXXXXXXXXXXX 
50XXXXXXXXXXXXXXXXX 



CHANNELS 10 

1 



1 2 3 4 5 

12345678901234567890123456789012345678901234567890 



CONTENTS 1. 1 1 1 1 11 1111114 3 1 

24671324925277104866432111123344475879301337560150 

28690015473255389571677832935004929789081291610120 

70652615438725664302854714164789679110584093597460 

66856629942783135483331165078129549753798806027850 



LOW-EDGE 



1. 1 

0998887766655444332221 100000 1 122233444556667788899 
06284062840628406284062840482604826048260482604826 



ENTRIES 
BIN WID 



26980 
.4000E-01 



ALL CHANNELS = .4134E+02 
MEAN VALUE = . 1525E+00 



UNDERFL 
R . M . 



done. 

krdriv: message: bye. 
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